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ABSTRACT 

Accurate mass determination of clusters of galaxies is crucial if they are to be used as cosmological probes. 
However, there are some discrepancies between cluster masses determined based on gravitational lensing, 
and X-ray observations assuming strict hydrostatic equilibrium (i.e., the equilibrium gas pressure is provided 
entirely by thermal pressure). Cosmological simulations suggest that turbulent gas motions remaining from 
hierarchical structure formation may provide a significant contribution to the equilibrium pressure in clusters. 
We analyze a sample of massive clusters of galaxies drawn from high resolution cosmological simulations, and 
find a significant contribution (20%— 45%) from non-thermal pressure near the center of relaxed clusters, and, 
in accord with previous studies, a minimum contribution at about 0.1 R w \ r , growing to about 3Q%-45% at the 
virial radius, R v [ r . Our results strongly suggest that relaxed clusters should have significant non-thermal support 
in their core region. As an example, we test the validity of strict hydrostatic equilibrium in the well-studied 
massive galaxy cluster Abell 1689 using the latest high resolution gravitational lensing and X-ray observations. 
We find a contribution of about 40% from non-thermal pressure within the core region of A 1689, suggesting 
an alternate explanation for the mass discrepancy: the strict hydrostatic equilibrium is not valid in this region. 

Subject headings: cosmology: theory-methods: numerical-gravitational lensing-X-rays: galaxies: clusters- 
galaxies: clusters: individual (A1689) 



1. INTRODUCTION 

Clusters of galaxies, the most massive virialized systems, 
form from the largest positive density fluctuations. The evo- 
lution of the abundance of these rare fluctuations are sensi- 
tive to the cosmological model. Also, the distribution of dark 
matter and gas in these systems provide a powerful test for 
our structure formation theories. Mass determinations based 
on X-ray observations customarily assume spherical symme- 
try and strict hydrostatic equilibrium, i.e., the gas pr essure is 
provi ded entirely by thermal pressure, = P t h (e.g.. lSarazinl 
1988). However, cosmological simulations suggest that even 
after equilibrium is established, a significant fraction of the 
pressure support against gravity comes from subsonic ran- 
dom gas m otion in clusters (Lau, Kravtsov & Nagai 2009; 
iMaieret all 120091: iFang et al.1120091: 1 Younger & Bryanll2007l: 
and references therein). 

Previous studies focused on how cluster mass deter- 
minations are influenced by non- thermal pressure support 
dZhang et alj l2010l; iM eneghetti et al| 120091: lLau et"aT] 120091; 
Lagana, de Souza & Kellerl 120091) . In this Letter, instead of 
determining the cluster mass, we focus on the dynamically 
important physical parameter of the intra-cluster gas (ICG): 
the pressure. We use a sample of massive clusters of galax- 
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ies drawn from high resolution cosmological simulations and 
quantify the contribution from non-thermal pressure in re- 
laxed clusters. 

As an example, we test the validity of the common as- 
sumption of strict hydrostatic equilibrium in Abell 1689. 
A 1689 is one of the most thoroughly investigated mas- 
sive galaxy clusters with a total mass of 1.5 x 10 15 
r'M w , located at a redshift of 0.183 (ICoeet al.1 120101: 
Peng et al.1 120091: iRiemer-Sejrensen et al.l 120091; iLemze et al 



2008; Umetsu & Broadhurst 2008; and references therein). It 



has been found, under the assumption of spherical symme- 
try and strict hydrostatic equilibrium, that in the central re- 
gion of A 1689 the mass derived from X-ray observations is 
significantly lower than that inferred from gravitational lens- 
m easurements dj| iemer-S0 rensen et al.l 120091: [Peng et al. 
lAndersson & Madeiskill2004t and references therein). 
We quantify the contribution from non-thermal pressure in 
A1689 using the latest CHANDRA, SUZAKU and gravitational 
lensing observations, and make use of our results for sim- 
ulated clusters to interpret the observations. In the rest of 
the paper we assume a concordance cold dark matter (CDM) 
model with Q m = 0.3, J7a = 0.7, and h = 0.7, where h is defined 
as Ho = 100 h km s _1 Mpc _I . Errors, error bars and dashed 
lines represent 1 a confidence levels unless otherwise stated. 

2. PRESSURE IN SIMULATED CLUSTERS 

We derive the hydrostatic and thermal gas pressure pro- 
files, P t h and P/ le , for massive relaxed clusters drawn from 
cosmological adaptive mesh refinement (AM R) simulations 
perfo rmed with the cosmological code ENZO (lO'Shea et al. I 
2004) assuming a spatially flat cosmological model very sim- 
ilar to the concordance model. The AMR simulations were 
adiabatic in the sense that no heating, cooling, or fee dback 
were included (for details see lYounger & Brvanl l2007). The 
highest resolution was about 25 kp c at r = 0.01 j? v jr (wher e 
the virial radius, R^, is defined as in lBrvan & Normanj l998), 
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FIG. 1 . — From top to bottom: radial profiles of dark matter and gas density 
(upper and lower curves) in units of critical density, p c ; gas temperature (in 
keV); hydrostatic and thermal gas pressure (upper and lower curves, in dyn 
cm -2 ); and pressure ratios, P t h/Phe, f° r massive simulated clusters. Solid 
lines and points with error bars represent best fit models and data points for 
our relaxed clusters (blue: CL1; green: CL2), and for a cluster with a non- 
relaxed core (red: CL3). We also show pressure ratios using gas pressure data 
points (blue plus signs, green stars and red crosses; see text for details). 

We selected relaxed clusters from the ten most massive clus- 
ters (1-2 xlO 15 M©) based on having a smooth spherically 
averaged density pro file and no sign of recent major merger 
events (for details see Molnar et al. 2009|). 

We show radial profiles for the two relaxed clusters (CL1 
and CL2) and one cluster with a non-relaxed core (but oth- 
erwise relaxed, CL3) in Figure Q] In this figure, the error 
bars represent the rms of the density and temperature varia- 
tions due to angular averaging. In the density and tempera- 
ture plots the solid lines represent best fit models to the AMR 
data points assuming double f3 models for the electron den- 
sity distribution, n e = «i/(l + r 1 / r\) l5Pi +n 2 /(l +r 2 /r|) L5A , 
and T e = T [l + T l exp{-(r/b l ) ai }]/(l + r/b 2 T 2 for the tem- 
perature. We fixed a 2 = 1 .6, which is an averag e value sug- 
gested by simulations (Youn ger & BrvarJ 12007). We used 
Pd oc l/[(r/di) l+a (l+r/di) 3 ~ a ] for the dark matter density 
profiles. 

We derive the hydrostatic gas pressure, P/ le , using the equa- 
tion of hydrostatic equilibrium assuming spherical symmetry, 



dPhAr) _ P g (r)GM(r) 
dr r 2 



(1) 



where p g {r) and M(r) are the gas density and the cumulative 
total mass within the 3-dimensional (3D) radius, r and G is 



the gravitational constant. We numerically integrate Equa- 
tion Q] using our fits. The accretion shock maintains a finite 
value for at R vlI , which we derive iteratively: we integrate 
Equation[T]inward from R^ demanding that the falloff of the 
pressure at large radii is a smooth function of the radius. Note 
that, since the pressure drops more than three orders of magni- 
tude from the cluster center to R v \ r , the pressure in the central 
region is insensitive to the choice of Pi, e at R v i r (within reason- 
able limits). 

The resulting hydrostatic and thermal gas pressure profiles 
are shown in Figure Q] In the 3rd panel the points show P,i,. 
The solid lines connecting the points are not fits but derived 
from the density and temperature fits (and serve as a consis- 
tency test). The sold lines above these lines represent the hy- 
drostatic equilibrium pressure. Note that, in the core region of 
CL3, the positive red error bars for the temperature are much 
larger than the negative ones. This is an indication of large 
positive deviations from the median in these radial bins, and 
a consequence of a non-relaxed core in CL3. The pressure ra- 
tio, R pr = P,h/Pi, e , profiles are also shown (bottom panel). In 
this panel, the errors (blue, green and red dashed lines) repre- 
sent the average deviations of pressure ratio points (plus signs, 
crosses and stars) from the mean (solid lines). We use these 
errors since they are larger than the statistical errors associ- 
ated with our smooth fitting functions (except at the center of 
the cluster with the non-relaxed core, CL3). 

3. TESTING HYDROSTATIC EQUILIBRIUM IN SIMULATED 
CLUSTERS 

We test the assumption of strict hydrostatic equilibrium in 
our simulated clusters using the pressure ratio profiles: R pr = 
Pih/Phe- These profiles in our relaxed simulated clusters show 
a similar trend (Figure[U bottom panel, blue and green lines). 
In each cluster, R pr is small near the center (r <, 0.05 /?„,), 
close to unity in 0.05 /? VJ > ^ r ^ 0.27?,,,y, and decreases in the 
outer regions (r ^ 0.2 /?,,,>). Even in our cluster with a non- 
relaxed core, CL3, the average R pl shows a similar trend (red 
lines). 

Quantitatively, we find that the contribution from non- 
thermal pressure, P„,/, = Pf, e -P t h, is decreasing from 20%- 
45% near the cluster center (0.01/? v ir), reaches a minimum of 
5%-30% at about 0.1/?vir, and increases up to about 30%- 
45% at R vlr . In the region of overlap with other simulations, 
r>0.1R V j r , similar contributions have been found in clusters 
from P„,h'. 5%-15% at about 0.1j? v i r , increasing with radius 
to 20%^10% at R vk (lLau et alJl2009l) . Our results for non- 
thermal pressure s uppor t are also consistent with those of 
Youn ger & Bryanl ((2007). Thus we find that in the central 
regions of our massive clusters, P,i, is significantly less than 
P he (Figure [B- 

In our adiabatic simulations, this non-thermal pressure sup- 
port is provided by subsonic random gas motion, i .e., Pd vn = 
P nt h, where Pd yn is the dynamical pressure (e.g., lLau et al] 
2009). We show the average dynamical pressure ratio, 
(Pdyn/Phe), profile in Figure [3] (black lines). We estimate the 
errors in (Pd yn /Phe) a&Min{(R pr —lcr)i} and Max{(R pr +la)i} 
(black dashed lines). We conclude that in the central regions 
of massive simulated clusters there is a significant contribu- 
tion from dynamical pressure. 

Based on our results, we find that there is no need for feed- 
back from a central AGN for the strict hydrostatic equilib- 
rium to break down in the central regions of clusters. This 
break down is due to subsonic random gas motion, which is a 
direct consequence of hierarchical structure formation. Near 
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FIG. 2. — Pressure ratio profiles, P t h/Phe derived for simulated massive 
clusters CL1, CL2 and CL3 (code is the same for the colored solid liens 
and points (as in Figure [T}, and for A1689 using the latest NFW lensing 
model (cyan lines, same as in Figuref5J. We also show the average dynamical 
pressure ratios, (Pj y „/Pi le ), for simulated relaxed clusters CL1 and CL2 with 
errors (black lines). See text for details and the definition of errors (dashed 
lines). 

R v i T , the contribution from P„,h may reach 3Q%-A5%. This is 
probably due to more recent and still ongoing slow accretion 
since the outer regions of clusters have not reached equilib- 
rium due to the large soun d crossing time, about 1 Gyr (for a 
more detailed analysis see Kawah arada et al.ll2010l) . Our re- 
sults for simulated clusters strongly suggest that even relaxed 
clusters should have significant non-thermal support in their 
core regions. 

4. TESTING HYDROSTATIC EQUILIBRIUM IN A1689 

As an example, we determine the pressure profiles in the 
well studied massive galaxy cluster A1689. We derive f,/, 
from our analysis of the three longest publicly available 
Chandra ACIS-I observations of A1689 (ObsIDs: 5004, 
6930 and 7289). After following standard ACIS data prepa- 
ration 9 , the total exposure time of 171 ksec was reduced to 
128.7 ksec. We fixed the redshift and the photoelectric ab- 
sorption at the galactic value, and used n = n e + nn, where 
n e and are the ele ctron and hydrog en densities assuming 
n e /riH = 1.1737 (as in lPeng e~ al. 2009). Since these observa- 
tions were taken in the Very FAINT telemetry mode, while the 
FAINT mode was used for the background observations, we 
chose to use local background. Assuming spherical symme- 
try, we determined the de-projected 3D densities and electron 
temperatures, T e , out to 0.6 for each shell fitting all shells 
and data sets simultaneously using XSPEC-12.5 applying C- 
statistic (Figure |3j black points with error bars). Our results 
are consistent with those Peng et al. (2009; see their Figure 
11). We calculate the thermal pressure in A1689 using the 
ideal gas law. For r > 0.6 R v i T we show the results for offse t 
1 from the SUZAKU observations (Kawaharada et al. 2010), 
which we use in our analysis. We also show the results for 
offset 2 (green points) to represent the similar results from 
offsets 2, 3 and 4 of SUZAKU. The lower resolution SUZAKU 
results agree with those from CHANDRA for r < 0.6 R vlr , thus 
we use only the last radial data points from SUZAKU. We use 
offset 1 because the other offsets would give us even higher 
non-thermal pressure contribution due to the very low temper- 
atures for offsets 2, 3, and 4 relative to offset 1. 

We derive the hydrostatic gas pressure, numerically in- 
tegrating Equation[T]assuming the same gas density and tem- 

9 http://cxc.harvard.edu/ciao/guides/acis_data.html 
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FIG. 3. — Top 3 panels: radial profiles of electron density and temperature, 
and thermal pressure in A1689 derived from CHANDRA observations (black 
points with error bars), and adopted from SUZAKU observations (blue points 
with error bars; offsetl; Kawaharada et al. 2010). On the op 2 panels we also 
show the best fit models (black solid lines), and, as a reference, the last point 
of offset 2 of SUZAKU (green stars with error bars). In the 3rd panel, the 
black solid line is derived from the density and temperature fits. In this panel, 
we also show hydrostatic equilibrium pressure, P/, f . derived from gravita- 
tional lensing using an NFW and a non-parametric model for the total mass 
distribution (cyan and magenta solid lines). Pressure ratios, R pr = P,i,/Pi w , 
are shown in the bottom panel (color code is the same as in the 3rd panel; 
horizontal black dashed line represents ratio of unity). 



perature models as in Section |3] We find the best fit pa- 
rameters for the density: n\ = 3.0 x 10~ 2 ±1.1 x 10~ 3 cm~ 3 , 
n 2 = 6.0 x 10" 3 ±2.7 x 10" 4 cm- 3 , (ri.ft) = (7.5 x 10" 2 ± 
3.8 x 10" 3 , 1.01 ± 6.0 x 10" 2 ), (r 2 ,A,) = (3.1 x 10^8.8 x 
10 _3 ,0.90±0.18x 10" 3 ); and for the temperature: T = 17. 1± 
3.6keV, 7] = -0.48 ±0.12 keV, (b h ai) = (0.40 ±0.15, 1.3 ± 
0.73), and {b 2l a 2 ) = (2.8 ± 2.3,2.3 ± 2.2), where all scale 
radii, r\,r2-,b\,b2, are in Mpc (Figure [3] black solid lines). 
The fitted parameters determining the large scale tempera- 
ture profile, (&2,«2X have large errors due to the large un- 
certainty in the SUZAKU temperature measurement near R w \ r . 
We derive the gas density using p g = n e n e nip, w here mp is 
the p roton mass (we adopt /i e = 1.146; see e.g., iPeng et aT] 
2009). We show the thermal gas pressure, P,i„ derived from 
the best fit density and temperature profiles in Figure [3] (3rd 
panel; black line). As a consistency check, we also derive 
the thermal pressure directly from the deprojected density 
and temperature points (points with error bars, same panel). 
We determine the cumulative mass function, M(r) using the 
best-fit spherical Navarro, Fre nk & White (NFW) model 
dNavarro. Frenk & Whitd[T997l) based on the latest high res- 
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olution mass modeling of A1689 by ICoe et al.l (|2010). As a 
comparison, we also use M(r) based on a non-parametric de- 
projection of the 2D mass distribution derived from a joint 
strong and weak lensing analysis of Hubble Space Telescope 
Advanced Cam era for Surveys (ACS) and S UBARU observa- 
tions of A 1689 (Umetsu & Broadhurst 2008) assuming spher- 
ical symmetry (see also|Kawaharada et al .1 120101) . Our results 
for Ph e are shown in Figure [3] (panel 3). The errors for the 
NFW model were estimated using Monte Carlo simulations 
assuming a Gaussian probability distribution for the parame- 
ters centered on their respective best fit values. The errors for 
the non-parametric model were estimated from Monte Carlo 
simulations based on the full covar iance matrix of the l ensing 
convergence profile (for details see lUmetsu et al. 2011J. 

The pressure ratios derived from observations of A1689, 
are shown in Figure [3] (bottom panel). These ratios, based 
on NFW and a non-parametric model for M(r), are consistent 
with each other (Figure|3j cyan and magenta lines). Similarly 
to our simulated clusters, we find a significant contribution 
from P„,i, in the core region and near /? v i r (Figure [3). The 
contribution from P nth is about 40% ±10% in A 1689 out to 
about 0.1/? vh, which is somewhat smaller than that in simu- 
lated clusters, although, at the center, it is within the errors 
due to cluster-to-cluster scatter. Also, from this figure we see 
that the maximum of the average ratio in simulated clusters is 
located somewhat closer to the center than the maximum in 
A1689 (Figure[3]l. 

However, we expect that in the core region of A 1689 non- 
gravitational processes are also important due to feedback 
from active galactic nuclei dMcCarthy et al..ll2009t and ref- 
erences therein). In general, contributions to the pressure 
in clusters may co me from turbulen ce, magnetic fields and 
cosmic rays (e.g., IVazza et ail 120091) . An attempt to sepa- 
rate these non-thermal pressure components based on their 
assumed functional form derived from numerical simulations 
has been carried out by Lagana et al. (2009). However, as 
they pointed out, this decomposition strongly depends on the 
assumed functional forms of the components. Unfortunately, 
at this point, we cannot determine the dominant contribution 
to the non-thermal pressure support in A 1689 from observa- 
tions directly. 

5. CONCLUSION 

We have analyzed massive clusters of galaxies drawn from 
high resolution cosmological simulations and quantified the 
non-thermal pressure support in relaxed clusters with high 
resolution from r ss 0.01/? v ir out to R w i r . We have found a sig- 
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nificant contribution from non-thermal pressure in simulated 
clusters due to subsonic random gas motion: 2Q%-45% at 
r w 0.01/?vii- and 3Q%-45% at R v i T having a minimum support 
of 5%-30% at r w 0.1,/?™ (FigureQ}. Our results strongly sug- 
gest that relaxed clusters should have significant non-thermal 
support in their core region, and that this non-thermal pressure 
support should be taken into account when analyzing clusters. 

As a test case, we have quantified the non-thermal pres- 
sure support in the well studied galaxy cluster A1689. Based 
on our results for the thermal gas pressure and the hydrostatic 
equilibrium pressure determined from X-ray and gravitational 
lensing observations of A 1689, assuming spherical symme- 
try, we have found a significant, 40±10%, contribution from 
non-thermal pressure within 0. 1 R vll . We conclude that the 
mass discrepancy in the central region of A 1689 can be ex- 
plained if we assume that the strict hydrostatic equilibrium is 
not valid in this region. We need to test the assumption of 
hydrostatic equilibrium in more clusters to find out how com- 
mon this large amount of non-thermal pressure contribution 
in their core region. 

While our spherical models for A1689 take into account 
support from non-thermal gas pressure as suggested by CDM 
models, they pr edict a high concentra ti on parameter, t he bi- 
axi al models of iMorandi et all ( 120101) . ICorless et ail (120091) 
and lSereno et al.l (1 2006b do not take into account non-thermal 
pressure but provide a concentratio n parameter con sistent 
with the predictions of CDM models. iPeng et all (120091) have 
found that a prolate gas distribution could solve the mass dis- 
crepancy, but it overestimates the total mass at large radii sig- 
nificantly, and implies a larger ellipticity than predicted by 
CDM models. Although all of these models solve the mass 
discrepancy in A 1689, neither of them is fully consistent with 
all predictions of CDM models. Our results suggest that a 
physical cluster model for A 1689 with a triaxial mass distri- 
bution including support from non-thermal pressure might be 
fully consistent with all observations and the predictions of 
CDM models. 
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